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ree PRODUCT ION 


A. BACKGROUND 

During the late 1800's and early 1900's, researchers 
memonteca finding periodic variations in mental, physical, 
and emotional capabilities of humans. These reported cycles 
form the basis of biorhythm theory. 

The proponents of Biorhythm theory have hypothesized 
miateethe mental cycle 1s 33 days long, the emotional cycle 
momzo days long, and the physical cycle is 23 days long. 
piGenermore they have suggested that these cycles begin on 
maemcday Of birth and continue throughout life with such 
Pmecision that an individual can calculate his position with 


megard to any cycle. 


Eee Se ORHYTHM DISCOVERY 

Thommen [1973] has described the work of early biorhythm 
researchers and the progression from discovery to present. 
The discovery of each cycle as presented by Thommen is out- 
ime De low. 

The initial development began about the turn of the 
century with the work of Dr. Herman Swoboda. His original 
investigations dealt with and culminated in his postulation 
of two basic cycles; one 23-day and one 28-day. Swoboda was 
principally a psychologist and was more interested in inves- 
PlifMemmoewiietier Man's feelings and actions were also in- 


fluenced by these cyclic phenomena. 





Gomeurrentty and independently, Dr. Wilhelm Fliess ob- 
served a 23 and 28 day cycle in the outbreak of fevers and 
deaths. Dr. Fliess had a well founded background in mathe- 
matics and statistics. He presented an in depth analysis 
@rehis work in a 564 page book entitled The Course of Life 
(Thommen, 1973]. 

Dieting ene 1920°s Dr, Alfred Teltsher, who had a doctor- 
mime in cngineering, analyzed the performance records of 
§,000 high school and college students. From his analysis 
he concluded that the high and low fluctuations in perform- 
mcm ollowed a definite 33-day cycle. From 1928 to 1932, 
tee Rextord Hersey and Dr. Michael Bennet conducted a study 
of emotions of workers in a railroad shop for the Pennsylvania 
Road. An analysis of this study revealed a cyclic variation 
mieemotlons with a period of 33 to 36 days [Thommen, 1973]. 

From the works of these early researchers the following 
descriptions of the three biorhythms have been postulated. 

It has been hypothesized that the physical cycle is 
fmeemty-three days long. The first 11% days of the cycle are 
mies pOsttiave half and the last 11% days are negative. The 
wie sewienre a transition from positive to negative or nega- 
Grve to positive CeGuremiomueriea a Critical day. The crit- 
here period uUsWally includes 12 hours on either side of the 
poet Or transition. Theoretically, during the positive half 
Semememecycle an individual is stronger and is better suited 
for heavy physical work and athletic contests. During the 


Seep eriod the individual is more subject to physical 





aziments. As such, it has been suggested that during the 
negative half individuals should avoid over exertion and 
@etempt to use this period for rest. [Willis, 1973]. 

The emotional cycle is 28-days long with each half cycle 
iMm—days bong. During the first haif (positive half) indi- 
fauals tend to be Emotionally more stable than the negative 
Zoe. «=Maxamum instability occurs during the critical period. 
[Willis, 1973]. 

The intellectual cycle is 33-days long and is divided 
into two 164 day segments. During the first period, the 
fme@uvrdual 1s mentally more alert than during the second 
half. Again the critical day is the period of minimum in- 


Meamkectual stability [Willis, 1973]. 


ee RECENT RESEARCH 

Posunpreecent researen in the area of biorhythms has been 
the correlation of accidents, deaths and illnesses with 
various phases of the biorhythmic cycles. Three of these 
eitics as reported by Willis [Willis, 1973] are outlined 
below. 

In Japan, the Ohmi Railroad Company records showed that 
pemoolmeaccidents during the past year 92 or 28% occurred on 
a critical day of the operator. If these accidents had 
Seeurred randomly, we could expect 20.4% to fall on a single, 
femobe, Or triple critical days. The probability of this 
Seetiememece OF a more extreme occurrence (greater than 92) is 
-001 assuming complete random accident occurrences. The Ohmi 


Railroad has instructed its operators in biorhythm theory and 





moarpnemtead a procedure that informs the operator of his 
Beoripethnmic state. Since adopting this procedure, the 
Negeaiama Service Center has logged over 2,000,000 accident 
free kilometers. [Willis, 1973] 

Miatysiseot patwent records at a medical facility as to 
Seats, accidents and medical incidents revealed that 50% 
ere the accidents and deaths occurred on a critical day. 

Bae percent of heart attacks and other such medical inci- 
dents occurred‘on a critical day. If these incidents were 
@astributed unaformly, only 20.4% would be expected to occur 
on critical days. [Willis, 1973] 

A survey of 100 single car fatalities in which the 
feeever was killed yielded the contents of Table I. For 
comparison the expected number in each category assuming 


complete random occurrences is also shown. 


Critical 

Days Obs # Expected | Race, ObS/ EXD. 
PEI 1 2035 SAE 

PE 6 LBs 0 we 

EI 4 .40 1Oz0 

PI Z ~49 4.08 

Pp 16 1.0 de Mt 

It ie() Sal 1.96 

EE 8 6.1 oa 


Pays heal, fo MoOt dona I , =e kee gill 
fete st. “Observed and Expected Accidents on Critical Days. 


DeeOsIECTIVES 
To date, most research in biorhythms has been in the 


Silat ysiseot iteld accident data. This type of analvsis seeks 
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to correlate individual accidents to various phases of the 
BrigGay time ctirve.  oince this line of investigation typically 
uses only one data point per individual, it is difficult to 
derive quantitative measures of variation within an indi- 
vidual. Another disadvantage to field research of this 
nature resides in the fact that the analysis of the data 
usually requires a subjective evaluation as to who or what 
mawecad tne incident. The present experiment seeks to over- 
come these disadvantages by measuring a performance parame- 
mer of an individual over a long period of time in a 
Semtrolled laboratory setting. Therefore, the primary ob- 
meen ve Of this research was to attempt to determine the 
presence or absence of pemrodiee rs UectUdatrens 1n performance 
on a laboratory task as biorhythm theory suggests. Secondly, 
ma periodic fluctuation was observed, to attempt to deter- 


mine the periods and amplitudes of such fluctuations. 
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Mire ood TAL METHOD 


A. EQUIPMENT 

The equipment consisted of a panel upon which the digits 
Mmecnrough $8 could be displayed and 8 response buttons. Each 
response button was related to one and only one of the digits 
1 through 8. This relationship was kept constant throughout 
Mmienwexperiment. Ihe ordering of the buttons was such that 
Powe One COrresponded to the digits 1 through 4 in ascending 
der trom left to right. Row two corresponded to the digits 
meerroupin Ss in ascending order left to right. The physicai 
relationship of the stimulus panel and response buttons is 


shown in Figure 1. 


Ste mmMuleus 
Display 


Response 
Buttons 





Frgume 1. 
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pe “OUBIECTS 

Three unpaid volunteer subjects were utilized in the 
Miecescut, experiment. All were male officer students at the 
Naval Postgraduate School. Subject one was 36 years old, 
subject two was 35, and subject three was 30. During the 
course of the experiment, the subjects were not given any 


biorhythmic data. 


oe TASK 

The experimental task required the subject to respond 
Pemoumckly as possible to the displayed stimuli. The re- 
Sponse required was depressing the button corresponding to 
the displayed digit. Only one response was allowed for any 
Ome Stimuli. The stimulus digits were displayed at a rate 
Gemone per second. Ihe ordering of the digits was random 
as determined by a random number table. This ordering was 
Changed daily and the ordering was not known by the subject. 
Each day, before data collection began, each subject was 


allowed a ten second warm up period. 


D. DATA COLLECTION 
The data collection was conducted at the same time each 
Baayveimeorder to nullify any circadian variations. The col- 
Bocr RON af data was carried out in three phases. During 
phase one, only the digits one and two were presented in a 
fandom sequence, this phase continued until five response 
Puiesmewenc recorded. During the second phase, the digits 


one through four were displayed in a random order. This was 


aS 





continued until five response times were recorded. During 
phase three, the digits one through eight were displayed in 
a random order and five response times were recorded. Be- 
fore the beginning of each phase, the subject was informed 
as to which digits would be displayed and that each digit 


was equally likely to appear during each display period. 
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III. DATA ANALYSIS 


A. COMPUTING MOVEMENT TIME AND INFORMATION PROCESSING RATE 

Experiments carried out by Fitts, Peterson and Wolpe 
moos | provided G€vidence which suggests that reaction time 
memaeemction of the probabilities of the stimuli. (i.e., 
ie more likely an event is to occur, less is the time re- 
marred tO react to it.} tTIhe functional form they observed 
was 


(1) RT get DE 


where 


RT = reaction time 


H = uncertainty associated with the stimuli (bits 
Sirona lO) 
b = seconds/bit of information processed 


1/b = information processing rate (bits per second) 
a = movement time. 
H is a term developed in communication theory by Shannon and 


Weaver [1949]. From their work, 


(2) Hi 


til 
n mo 
_— 


P (i) log, p(i) bits/element. 


th 


P(i) PRODdbmii tyeot the i Stimuli occurring. 
For n equally likely stimuli (eg(2)) becomes 
(3) He log, We 
For present work, three tasks were administered sequen- 


tially with five repetitions of each task. The task involved 


15 


making a known response to an rece can SCimelen. = ine wransit 
task utilized two equally likely stimuli (one bit), the 
second utilized four equally likely stimuli (two bits) and 
the third used eight equally likely stimuli (three bits). 
The fifteen reaction times were then used in a linear 
meeression model to produce a least square error estimate 


of a and b. The equations used were: 


on 


5 15 
- ae 
1 7? i=l 


= 


i, 
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r feo 
mee i Me 
rm nie 

bk 
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1 
where 
ees Dit S Or onl Ormlat ton 
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ze 6 
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number of data points. 


15 lS 


(5) oo - b 


RB. DETERMINING SIGNIFICANT VARIATION 

Pemelemecotmelistion Of the data collection phase, data for 
each subject consisted of three sets of measurements. Set 
Ole was tne average reaction time. Each point in this set 


was the average of the fifteen reaction times recorded during 
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a daily session. The distributional form for the two bit 
reaction time was found to be approximately Normal (u,o7). 
The assumption was made that on any day the distribution of 
the one bit decision was Normal (yu - c,o*) and the distri- 
bution of the three bit decision was Normal (u + c,o*) where 
c is a constant. The constant c was estimated as the inverse 
of the information processing rate. The constant c was added 
to the one bit data and subtracted from the three bit data 

to yield fifteen data points which were distributed Normal 
(u,07). The variance was estimated using the unbiased 


estimator, 


2 1 a= 2 
fe) ~ sai 2 ee) 5 
n=1 
where 
xX = mean 
n = number of data points 
x, = value of nee data point. 


The mean of the distribution was estimated by 


igs) 
a, oe 
=] nN 


x = 


=O Noe 


n 


With the estimated value of o*, the confidence interval for 
@e= .05 was calculated using the Student t distribution. 
The relation used was L = 2 ae (1-a/2)(S?/n), where 

n - 1 = degrees of freedom 


number of data points 


nN 


(1-a/2) 


APinertetOn to distribution. 


Ly 





ieee, nemaca the tollowing confidence intervals centered 


around x. 
Subject one: L = ,036. 
Subject two: L = .052. 


Subject three: L = .044, 

The grand mean (X_) for each subject was computed and 
the data was analyzed to see if any data points lie outside 
the interval X, = +L/2. In all cases this was found to be 
SO, indicating a significant variation in daily performance. 
This lead to the Fourier Analysis of the data set described 
im the next section. It was also decided to conduct a 
Fourier Analysis of the movement time (MT) data and the in- 


formation processing rate (IPR) data. 


C. FREQUENCY ANALYSIS OF DATA 

In order to determine periodicity within the data set 
it was decided to use a Fast Fourier Transform (FFT). Since 
this Sperictronescaquired Equal interval data, this necessi- 
tated generating data points for days when subjects were not 
available. This was done by a linear interpolation between 
the point immediately preceding and the point immediately 
following the interval with no data. This procedure was 
necessary for all weekend days and holidays. Although the 
exact effect of the manicien of data is not known, it was 
felt that most effects would be restricted to frequencies 
Peuieperi10ds Of 5 £0 7 days. The FFI indicated periodicities 


Within each data set. On this basis, the function described 
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by the data set was approximated by a Fourier Series of the 


form: 
a (NN) P -_ TL(NN)P 
Y(NN) = —+ f (a_ cos + bo sin ——— ), 
2 a n N n N 
n=] 
where 
Ne=3s (number of data points) 
NNee O,1,2, <-- {2N-1 (days) 
P= harmonic number 
Y(NN) = value of function on day NN. 


An existing subroutine from the Naval Postgraduate 
School computer library was used to compute the coefficients 
a. and bo The subroutine is attached as Appendix M. The 
root mean square value of a, and DL was computed to give 
mre amplitude cf zee [MniGimern cl crcacen tO get a bettcr 
mepEoximation of the frequency envelope at the lower fre- 
Gmencies, the data Set was reduced in increments of one data 
point until the final data set had only 46 points. This 
operation is valid under the assumption that the amplitude 
of the various frequency components do not change over time. 
The frequency envelopes of reaction time, movement time, and 
information processing rate are shown in Figures 2 through 
nO 

The phase of each harmonic can be found by using the 


following relation: 


a 
n 

6 = arctan Bo? 
n 


where a, and bo are the Fourier Coefficients. 
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The phase angle @ for periods with significant amplitudes 


are shown in Table II. 


D. LEVEL OF SIGNIFICANCE IN FREQUENCY DATA 
Under the Null Hypothesis, that is that f(t) is a Gaussian 
Process, the coefficients of the sine and cosine functions in 


the Fourier series are assumed to be independent, identically 


distributed Normal with mean of 0 and variance of o%. Since 
2 yay 2 
AL an : or 
= . th : 
A, ="amplitude of n harmonic 
then 
nee il 
AL exp Gaze 
o* can be estimated by the unbiased estimator 
2 1 s 2 
ave oe 2 en 
n=l 
where 
N = number of observations of a 


~ 
Since WW is assumed to be zero, this reduces to 
4 2 

Ned 2. “fn - 

nel 
The values of the S* computed are shown in Table II. 

inewpmobabilaty (a) an amplitude greater than or equal 
to A, would arise from a Gaussian Process is found from the 
relation 


= Ne aS 
a = © e . 


Zo 





In the original set of 70 points, there are 6 harmonics 
Prttecoula Dewsienaticant in the» range of 10 - 35 days. The 
confidence level that can be associated with the amplitude 


of one of these harmonics is: 


=e J 25° 
a=l1-(l-e 


These values of a for the peaks in the frequency envelopes 


are shown in Table II. 
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IV. RESULTS AND CONCLUSIONS 


Table II gives the data derived for the most obvious 
peaks in the various frequency envelopes. Of the 14 periods 
with amplitudes significant at the .1 level, nine are with- 
in one day of half of one of the three biorhythm cycles. 
meme robabitity Of this occurrence or a more extreme occur- 
rence is .015, assuming a uniform distribution of signifi- 
cant amplitudes. 

The significance of periods with one half the biorhythmic 
periods can be seen by examining the critical - non critical 
day cycle. There are two critical days within each period. 
Hence, the period of a critical day - non critical day cycle 
meone half the biorhythmic period. 

Although all subjects have significant periods within 
one day of 11.5, 14, and 16.5, an examination of the phases 
associated with each indicated that the periodic phenomena 
apparently arise from different sources. Since all subjects 
performed the same task each day, it can be assumed that the 
cyclic phenomena was not induced by the experiment. 

During the analysis of the data, it was found that two 
unknown effects were present. The first was the interpola- 
tion of missing data points. This effect can be eliminated 
PpetistiniaoncOntinuous availability of subjects. The other 
effect was the individual's environment. Outside of the 
test situation, each individual had a unique environment 


which may have affected his performance. It is suggested 


Sl 





Saat im tUrther investigations an attempt be made to obtain 
pubpyects With a Uniform environment during the course of the 
Seoeriment. 

An experiment involving a similar procedure but with 
subjects drawn from a population where extraneous influences 
can be controlled (e.g. prisons, Naval Academy) may promote 
obtaining information regarding the true influence of bio- 
rhythms. Also a larger population of subjects would allow 
estimates to be made of the variation in the influence of 


meonrnythms in the general population. 
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APPENDIX A 


Computation of probability of a critical day assuming 
independence and uniform distribution of critical days. 

1. Since two critical days occur during each cycle, if 
a large number (N) of people are sampled on any given day, 
there should be (2/23) (N) in a critical physical day; 
mee7oy(N) in a critical emotional day; and (2/33)(N) ina 
ermctical mental day. 


In terms of probabilities: 


-087 (probability of a physical critical day) 


Ee 

Pp 

ee = .071 (probability of an emotional critical day) 

Pen evoleprobaotlaty Of an intellectual critical day). 


Pe Ccomucing Ene number of Single critical days. 


(a) Physical - Et Bott ED = .076 
fo Emotional = ee = .061 
(c) Mental = = .051. 


he ee 
Total probability of single critical day = .188. 


3. Computing the number of double critical days. 


(a) PG&E = Eo) Oe Ose) = .0058 
(b) P & M = (P) (PGP) = .0049 
(c) M& E = 


ea) ea Ge) = ,0040. 
fotaleprotaprlity of double critical day = .0147. 
4. Computing the probability of a triple critical day 


Prriple = (Pe) (Py) (Pp) = -00038. 
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Weaetorvieprondamitity Of having a critical day (single, 


double, or triple) 


P single R “deuble ° Periple = oS 
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SUBJECT ONE 


Average Reaction Time 
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SUBJECT ONE 


Movement Time 
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One Reaction Time 
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SUBJECT TWO 


Movement Time 


“nteenmpolated Data Points 
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SUBJECT TWO 


Processing Kate 
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APPENDIX H 


SUBJECT THREE 


Date Average Reaction Time Date Average Reaction Time 
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APPENDIX K 


SUBJECT ONE - BIORHYTHMIC DATA 
Bain Dajieee2/ 10/38 


Number of daysinto each cycle on indicated date. 


Date baby sta I Emotional Mental Date Physical Emotional Mental 


4-30 4.0 ILS ae ae) 1) 6-4 or 2070 Saat 
s- I ay. 14.0 5020 G5 aber) Ms MU Sa0 
a-Z 6.0 lhsy.0 Siw G- 6 hess 0 22.0 S50 
a- 3 120) ILS SZ G= 7 We) 2 oe0 ie) 
5-4 Sa 0 SS) GizG Zw) 24.0 2.0 
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5-6 ko. Ie) we Za) Gi 10 Yao «Ni 260 4.0 
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Sa) Loe PETE Alt a) Go 13 ZU) 1.0 (acl 
ae) «=i 4.0 ZW) G0 6-14 S70) Zit &.0 
sel §€615.0 24.0 Pv) 6- 15 4.0 5750 2h 
fie 16.0 Za 8.0 O=A'G 9.0 4.0 Orso) 
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aoe 8621.0 Lal) Ss Ge Zall IL 2b 13.20 
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Says 0 Sca0) (EAMG dbs Ge 26 Sef 14.0 200 
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APPENDIX L 


SUBJECT TWO - BIORHYTHMIC DATA 


Birth Date: 8/15/39 1 i - 
Bees / TGA. ¢ 
Number of days into each cycle on indicated date. 


Date Physical Emotional Mental Date Physical Emotional Mental 
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> - 2) i fae) ASU Zion 0 6-24 WS). ZU 0 27.0 
ja ars | 14.0 Zio U Oe 2 14.0 2 Age 28.0 
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APPENDIX M 


Supeeel LMREB - BIORHYTHMIC DATA 


Bipth Date: 1/15/44 
Number of days into each cycle on indicated date. 


Date Physical Emotional Mental Date Biisieal Emotional Mental 
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29 By (0) OZ O20 poe Lay yee 8) Ee, 
Se 0 4.0 LAO Soe Lied hoy at oO rer, 
>= epi 220 7.0 0 (ia 170 Ie aU 
Ge L 6.0 >, 0 el) i= 0 1350 Zor 0 4.0 
OZ 12a) 14.0 Ov (a7 oso JAAN 30 52.0) 
Ge.5 ou ae 0 4.0 aS ZrO 2250 6.0 
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FORTRAN IV C LEVEL 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
ee, 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
003i 
0032 
O03 3 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 


0048. 


13 


23 


251 
ee 


2$2 


APPENDIX N 


Ze FOURIA DATE = 


SUBROUTINE FOURIA(LP, Y,M,AHC,SYS,A,B,E) 
DIMENSION Y(1) ,A(1),E(1) ,E(1) 

L=LP 

LM=L-1 

IL = 2 

N = L/2 

NFI=N=L 

IF(MCD(L,2).NE.0) IL=1 

ME=M 

IF(M.LE.0.CR.M.GT.N) ME=N 

SUMYS = 0 

Hor ta. 1 
SUMYS = 
SYS = § 
é 
C 
S 
FL = L 
FN = FL/2. 

CCN1 = 6.283185/FL 
C1=CCS(CCN1) 
S1=SIN(CCN1) 

CCN2 = 2./FL 


74211 


= 1,L 
SUMYS + Y(I) ** 2 
UMYS 


One il 


IM = 1 
D) 28 J = 1,NP1 

e= J = 1 

IF (J.EQ.NF1) IM-IM+tIL 
Bie = 0. 

Ul = 0. 

LN = L 

DO 23 K = 1,LM1 

UO = Y(LN) + Q*U1 - U2 

U2 = Ul 

U1 = UO 

LN = LN - 1 

AAA = CCN2 * (Y(1) + C* Ul - U2) 
cer Gh-0) GO TO 251 
AHO = AAA/2. 

SUMC = AHC * AAA 

GO TO 26 

eae. 3) GO TO 252 
A(N) = AAA/2. 

B(N) = 0. 

SUMC = SUMC + A(N)*AAA 

GO TO 254 

A(IF) = AAA 

B(IP) = CCN2 * § * Ul 
SUMC = SUMC + A(IP) ** 2 + B(IP) ** 2 
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0049 254 E(IP) = SUMYS - FN * SUMC 
G 


0050 
0051 
0052 
0053 
0054 
0055 


0056 
0057 


C 
C 


26 


28 


IF(IP.EG.ME) GO TO 33 
ee= Gis G - S1 * S 

ee ei S + SI *.C 

Ge= GC 

G = C+G 

CONTINUE 


33 RETURN 
END 
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